Lack of cooling flow clusters at z > 0.5 
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The goal of this work is to study the incidence rate of "cooling flows" in the 
high redshift clusters using Chandra observations of z > 0.5 objects from a 
new large, X-ray selected catalog [1]. We find that only a very small fraction of 
high-z objects have cuspy X-ray brightness profiles, which is a characteristic 
feature of the cooling flow clusters at 2; 0. The observed lack of cooling flows 
is most likely a consequence of a higher rate of major mergers at z > 0.5. 

1 Introduction 

The central regions in a large fraction of low-redshift clusters are clearly 
affected by radiative cooling [2]. Some estimates put the fraction of such 
cooling flow clusters to > 70% (e.g., [3]). A recent by Bauer et al. suggests 
that the cooling flow fraction remains high to z '-^ 0.4 [4]. However, this 
work is based on the ROSAT All-Sky Survey cluster sample, and so it can 
be strongly affected by Malmquist bias (strongly over-luminous clusters are 
preferentially selected because of the high flux threshold). 

Any evolution in the cooling flow fraction, if detected, must be taken 
into account in detailed physical models of this phenomenon. We address 
this important question using a new distant cluster sample, derived from a 
sensitive survey based on the ROSAT pointed observations [1]. All objects 
were observed with Chandra, providing a uniform dataset which should be 
much less affected by selection effects than the previous samples. 

2 Definition of the "cooling fiow" cluster 

First of all, we need to choose a definition of the "cooling flow" cluster that 
can be efficiently applied to the X-ray data of various statistical quality. The 
most common definition is based on the estimated central cooling time: cool- 
ing flow clusters have icooi ^ tn (e.g., [3]). One could also use the mass depo- 
sition rate given by the standard cooling flow model [2] ; cooling flow clusters 
have M > (10 — 100) Mq yr~^ [3]. These definitions rely on spatially-resolved 
spectroscopic measurements which is a serious disadvantage for application 
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Fig. 1. X-ray surface brightness profiles typical for non-cooling flow {left; A40i) 
and cooling flow {right; A85) clusters. Solid lines show the model X-ray brightness 
corresponding to the best-fit gas density model (see [5] for details). 



at high redsliifts. The data of sufRcient quality to measure T{r) are available 
only for a small number of high-z objects. We, therefore, look for a definition 
based solely on the X-imaging data. 

At low redshifts, there is a clear connection between the presence of the 
cooling flow and the X-ray morphology. Clusters with tcooi ^ tn have X-ray 
brightness profiles with flat cores while those with icooi ^ tn have character- 
istic central cusps in the X-ray brightness distribution (Fig. 1). The central 
cusp can be characterized by the power-law index of the gas density profile, 
a = d\ogpg/d\ogr. For uniformity, the radius at which a is measured should 
be chosen at a fixed fraction of the cluster virial radius. This radius should 
be sufficiently small so that the effects of cooling are strong. At very small 
radii, however, the density profiles even in clusters with strong cooling fiows 
can flatten because of the outflows from the central AGN (sec many papers 
in these proceedings). Empirically, a good choice is r = 0.04i?5oo,^ and so we 
define the "cuspiness parameter" , a, as 

^^dlogp, r = 0.04i?500 (1) 

alogr 

Cuspiness can be measured by fitting a realistic 3-dimensional gas density 
model to the observed X-ray surface brightness (our procedure is described 
in [5]). Examples of the best-fit models are shown by solid fines in Fig. 1. Such 
modeling is feasible with moderate-exposure Chandra observations of high- 
redshift clusters. i?5oo can be estimated using low-scatter cluster mass proxies 
such as the average temperature (excluding the central cooling region). We 
use an even better proxy, the recently proposed Yx parameter [6], which is 



^i?50o is the radius at which the mean enclosed total mass overdensity is 500 
relative to the critical density at the object redshift. iisoo ~ 0.5Rv 



Lvir . 
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Fig. 2. The distribution of the cuspiness parameter in the low-z cluster sample. 
Arrows indicate the values for some well-known clusters. The boundary value of 
a = 0.5 approximately separates cooling flow and non-cooling flow clusters. 



remarkably insensitive to the cluster dynamical state and easily measured 
even in high-redshift clusters. 

Our low-redshift cluster sample is a flux-limited subsample of 48 objects 
from the HIFLUGCS catalog [7] , all with the archival Chandra observations. 
The distribution of the cuspiness parameter for these objects is shown in 
Fig. 2. Clearly, the value of a is closely connected to the more common cooling 
flow definitions. Clusters with a > 0.7 (e.g., A2065, A478, A2029, A2597, 
2A 0035, A133) are known to host strong cooling flows. The objects with 
Q. < 0.5 (e.g., A2163, A399, A119, A2256, A754) are famous non cooling 
flow clusters. The clusters in the range 0.5 < a < 0.7 (e.g., A2589, A3571) 
host weak cooling flows. Therefore, the cooling flow clusters are those with 
a > 0.5. Approximately 2/3 of the low-redshift sample (31 of 48 objects) 
have cuspiness above this value, in line with the previous estimates of the 
cooling flow incidence rate (e.g. [3]). 



3 High-redshift cluster sample 

Our high-redshift sample is derived from the recently completed 400 deg^ 
ROSAT PSPC survey (400d; [1]). This is the largest-area survey based on 
the ROSAT pointed observations. Clusters are detected as extended X-ray 
sources in the central 17.5' of the PSPC FOV and required to have fluxes 
fx > 1.4 X 10""'^^ erg s~^ cm~^. The X-ray sample is fully identified. It 
includes 266 optically confirmed clusters (95% of the X-ray candidate list). 
Spectroscopic redshifts are available for all objects. 




Fig. 3. (a) — Detection efficiency of tlie 400d survey for idealized /9-model clusters 
as a function of total flux and core-radius. Dotted line shows the flux limit of the 
400d catalog. Detection efficiency is reduced for Vc < 8" because the such clusters 
are hard do distinguish from the point sources. The efficiency is also small for 
objects with Tc > 150" because they are "lost" in the cosmic X-ray background. 
(b) Detection efficiency as a function of core-radius for = 2x 10~^^ erg s"'^ cm~^. 
Shaded histogram shows the distribution of core-radii in a low-redshift sample [8] 
scaled to 2; = 0.5. 



A subsample of the high-z 400d clusters has been observed with Chandra. 
The exposure times were chosen to yield at least 2000 source counts, which is 
sufficient to measure the average cluster temperature with a 15% uncertainty 
and accurately trace the surface brightness profile to r ~ Rboo- Chandra^ s 
angular resolution corresponds to a linear scale of < 8 kpc out to z = 1, 
fully sufficient to measure the cuspiness parameter. In this work, we use 400d 
clusters at z > 0.5, 20 in total. The typical mass of these objects corresponds 
to today's 4 keV clusters. 

The basic characteristics of the X-ray selection in the 400d survey have 
been extensively calibrated through exhaustive Monte-Carlo simulations 
(see [1] for details). The aspect most relevant for the present study is the 
sensitivity of the detection efficiency to the cluster size and structure. A 
precise two-dimensional map of the detection efficiency as a function of clus- 
ter size and core radius was derived for idealized /3-model clusters (Fig. 3a). 
The detection efficiency drops significantly only for Tc < 8" and Tc > 120" 
(see Fig. 3b which shows the slice through the detection probability map at 
fx — 2 X 10~^^ erg s~^ cm~^, just above the survey flux limit). The angular 
size range in which the 400d X-ray detection is sensitive encompasses the 
entire range of core-radii expected for the high-redshift clusters (c.f. shaded 
histogram in Fig. 3b). Therefore, the 400d selection will not bias the distri- 
bution of core-radii for /3-model clusters. 

The sensitivity of the 400d X-ray detection algorithm to cooling flow 
clusters with the cuspy X-ray brightness profiles requires a separate study. 
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Fig. 4. (a) Chandra images of the 400d clusters with z > 0.5. Note a high fraction 
of objects that show clear signs of a major merger, (b) Distribution of the cuspiness 
parameter for the 2 > 0.5 and 2 ~ samples. 



This issue was addressed by a separate set of Monte-Carlo simulations in 
which instead of /3-models, we used the real X-ray images of a complete 
sample of low-z clusters, scaled to different redshifts in the range 0.35 < z < 
0.8 (see [1] for details). A short summary of the results from these simulations 
is that there is no significant difference in the detection efficiency for the 
/3- model and cooling flow clusters (see, e.g., Fig. 16 in [1]). For example, 
Hydra-A {a = 1.24) at z = 0.45'' is detected with the probability 0.67; 
2A 0335 (a = 1.26) at z = 0.45 has pdot = 0.54; A2029 (a = 0.90) at z = 0.8 
has pdct — 0.69. These values are near the maximum efficiency for /3-model 
clusters of similar flux (Fig. 3b). Therefore, there should be no discrimination 
in the 400d survey against the objects similar to today's cooling flow clusters. 



4 Observed morphologies and cuspiness parameters of 
the high-redshift clusters 

Chandra images of the z > 0.5 objects from the 400d sample show a clear 
evolution of the cluster X-ray morphologies — at least 15 of 20 objects shows 
signs of an on-going major merger (Fig. 4a); the corresponding fraction in the 
low-redshift sample is < 30%. The same effect is apparent in the distribution 



''Redshifts here are chosen so that the observed fluxes would correspond to that 
in Fig. 3b, 2 x 10"" erg s"^ cm"^. 
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of the cuspiness parameter shown in Fig. 4b. Only 3 of 20 high- 2; clusters 
have a > 0.5 (i.e., above the boundary between cooling flow and non cooling 
clusters, see § 2), while in the low-z sample this fraction is 31 of 48 (65%). 
The are no clusters with a > 0.7 (strong cooling flows) in the z > 0.5 sample, 
while the fraction of such clusters at z ^ is 46% (22 of 48 objects). The 
statistical significance of the difference in the distribution corresponds to a 
random fluctuation probability of only P ~ 5 x 10~^. 

Our results provide a tantalizing evidence for a strong evolution in the inci- 
dence rate of the cluster cooling flows at z > 0.5. This is apparently related 
to the higher cluster merging rate, indeed expected at these rcdshifts in the 
Dark Energy dominated, cold dark matter cosmological models (e.g., [9]). 
The cluster cooling flows thus appear to be a relatively recent phenomenon, 
which becomes common only in the past 1/3 of the Hubble time. 
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